Inhomogeneous spin diffusion in traps with cold atoms 
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The spin diffusion and damped oscillations are studied in the collision of two spin polarized 
clouds of cold atoms with resonant interactions. The strong density dependence of the diffusion 
coefficient leads to inhomogeneous spin diffusion that changes from central to surface spin flow as 
the temperature increases. The inhomogeneity and the smaller finite trap size significantly reduce 
the spin diffusion rate at low temperatures. The resulting spin diffusion rates and spin drag at longer 
time scales are compatible with measurements at low to high temperatures for resonant attractive 
interactions but are incompatible with a metastable ferromagnetic phase. This does not exclude 
that the colliding clouds can evolve into a repulsive initial state which subsequently decays during 
the bounce and the initial damped oscillations. 
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Measurements of spin diffusion in resonantly interact- 
ing trapped atomic gases [J provide new understanding 
of the crossover physics and ferromagnetic (FM) phases. 
The transitions observed in radii and expansion energies 
versus repulsive interaction strength [2] were interpreted 
as a FM Stoner transition [3j although disputed in Ref. 
4\. A FM transition is predicted in several calculations 
5 with qualitatively similar transitions in observables 
but at a lower interaction strength as in the experiments. 
Recent experiments on colliding spin polarized gases and 
subsequent spin diffusion seem to exclude a metastable 
FM phase for resonant interactions Since the spin dif- 
fusion in traps occurs on longer time scales than the ini- 
tial bounce of the spin separated clouds and subsequent 
damped oscillations, a FM phase initially as claimed by 
Taylor et al. [l(| is not necessarily excluded. 

Resonant interactions are only characterized by an in- 
finite scattering length which can either be in the unitar- 
ity limit of strongly attractive interactions (BCS-BEC 
crossover) or in the repulsive (FM phase). These are 
also referred to as the lower and upper branches respec- 
tively of the multivalued scattering length (see e.g. Q). 
The upper branch is metastable and decays to the lower 
branch due to three-body interactions [11) . 

The diffusion coefficients measured in the traps seems 
to exclude a FM [l| yet the diffusion rate exceeds the cal- 
culated values at central densities [IH by up to an order 
of magnitude in particular at higher temperatures. At 
temperatures well above the Fermi temperature T^iTp 
the difference can be explained as a density inhomogene- 
ity effect jlH arising because the spin diffusion is faster 
at the lower surface densities and as a result the spin 
diffusion circulates faster avoiding the nondiffusive core. 
The spin diffusion calculations will here be extended to 
low temperatures T^JTp including the effects of density 
inhomogeneities. Comparison to the measured diffusion 
coefficients and the initial damped oscillations will then 
indicate whether a FM phase was present for resonant 
interactions. 



Defining the spin density n s — — ni and current 
j s = jf — jj,, the spin diffusion coefficient is defined 
from Fick's law. For inhomogeneous systems the diffu- 
sion equation leads to a modified form for Fick's law [l3j 
js = --DxV(n„/x), where x = dn s /d(fi^ - m) is the 
spin susceptibility. Inserting the equation of continuity 
for spins d t n s + Vj s = 0, we obtain the diffusion equation 
for the spin diffusive mode n s (r, t) — e^ t ^ T n s (r) 



V(D X V(n s / X )) = 0. 



(1) 



The axial spin diffusion rate r^ 1 was calculated in Ref. 
[l3| in the unitarity limit and at high temperatures, 
where D\ is independent of total density n = + 
At low temperatures D\ oc n 5 / 3 in the unitarit y li mit 
whereas D\ oc n in the dilute limit according to |l2| : in 
both cases \ oc n 1 ^ 3 . The various density dependences 
are conveniently parametrized by D\ = DqXqti p /tIq, 
where p is a general power dependence decreasing from 
p = 5/3 at low temperatures to p = at high tem- 
peratures. Do, xo and uq are the values in the trap 
center. At low temperatures the atomic cloud density 
n(r) = rio(l — r\/Ft?^_ — z 2 /it^) 3 / 2 , has a sharp surface at 
the Thomas- Fermi radii Rf — 2£f?/mu 1 Vl + /?. At the 
BCS-BEC crossover the universal constant is j3 ~ —0.54 
for unpolarized gases. 

The diffusion equation (flj is precisely the one analyzed 
for collective modes in the case p = 1, where analytical 
results exist The diffusive mode along the z-axis 

corresponds to the axial spin dipole mode with n s /x oc 
z, which by insertion in the diffusion equation (fTJ) gives 
t _1 = 'SDq/R^ independent of trap deformation A = 
u) z /oj_l. 

The diffusion equation (fTJ) results from using the vari- 
ational principle on 



Jd 3 rD X (W(n s /x)) 2 
Jd 3 rn 2 Jx 



(2) 



' heiselberg@mil.dk 



Varying the spin density n s on the right hand side gives 
an upper estimate for the spin diffusion rate, and its mini- 
mum determines both the spin density and diffusion rate. 
It is experienced in transport calculations [15] that varia- 
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tional results are surprisingly accurate even with the sim- 
plest ansatz that obeys conservation laws and symmetry. 
We shall also find this in our case with n s /x oc z (see 
Fig. (1)). We emphasize that the diffusion equation ([T]) 
and the variational equation @ differ in our low temper- 
ature case from Ref. [13| because D\ is not independent 
of density and also by the sharp density cutoff at the 
surface of the atomic cloud as will be explained shortly. 
Also the qualitative and quantitative change in diffusion 
is accurately captured in the variational calculations. 

A variational calculation of Eq. ([2]) with the more gen- 
eral form n s /x oc z(l + c z z 2 /R 2 . + c±r 2 L /R 2 L ) results in 
a vanishing curvature c± = in the radial direction for 
prolate traps, A -C 1. Generally we find that any such 
radial curvature increase the diffusion rate proportional 
to / A 2 due to the gradient in Eq. © , and therefore it 
vanishes for finite prolate traps. The axial curvature c z 
is also independent of trap deformation when A <C 1 , but 
does not vanish except at p — 1, where it changes sign 
(see Fig. 1). A full variational calculation with z A and 
higher order contributions in n s / z\ shows that these cor- 
rections are small and reduce the spin diffusion rate only 
slightly for p^,0. This is towards the high temperature 
limit where the cloud extends outside the Thomas-Fermi 
radii anyway as discussed below. 

The variational minimum yields the diffusion rate 
T var — $Do/R~)Ip, which we write in terms on a di- 
mensionless rate I p that contains the inhomogeneity of 
the trap diffusion. It generally depends on p except as 
mentioned when p — 1 where I p = 1 but is indepen- 
dent of trap deformation when A <C 1. For prolate traps 
we find that the ansatz n s cx z\ is a good approxima- 
tion for positive p as shown in Fig. 1, since the flow 
predominantly takes place along the symmetry axis near 
the center. Inserting this ansatz in Eq. @, we get the 
(upper) estimate for the rate or in terms of the dimen- 
sionless "inhomogeneity" factor 

j _3_ 5r(4)r(3p/2 + 2) 

v 3D t (3p + 2)r(7/2)r(3p/2 + 5/2) ' [) 

given in terms of the r-function. The importance of in- 
homogeneous spin diffusion is evident as shown in Fig. 1. 
The rate decreases with increasing power p because the 
density and therefore the spin diffusion decrease toward 
the surface of the trap. Comparing temperatures T ~ Tp 
where p ~ to low temperatures where p = 5/3, we find 
that the inhomogeneity factor is ~ 6 times smaller at low 
temperatures. 

The diffusion rate at low temperatures is qualitatively 
different from the high temperature case studied in Ref. 
[l3j . At low temperatures the densities vanish outside 
the Thomas-Fermi radii and therefore spin density and 
susceptibility also vanish. Consequently, the length scale 
in the radial direction is limited by R± . This is contrary 
to high temperatures |13| . where the gas extends over a 
length scale k = y^2T/mujf but the spin diffusion ex- 
tends over a much longer radial length scale ~ l^/^/X 
because the diffusion coefficient increase with decreasing 
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Figure 1. (Color online) Spin diffusion rate versus the power 
of the density dependence for prolate traps, A < 1 , The 
inhomogeneity factor I p = R^/3Dor is calculated for n s /x oc 
z (dashed curve) and variationally with n s /x = z(l+c z z 2 / R^) 
(full curve), where c z (dash-dotted curve) has been varied to 
minimize the rate of Eq. (f2]). 

density. Therefore diffusion dominates in the surface and 
tail of the distribution such that the spin flow circulates 
around the trap center. As result the diffusion rate is in- 
creased by a factor A -2 as compared to the low temper- 
ature result of Eq. ([3]). However, because the diffusion 
eventually becomes collisionless at very low densities in 
the tail, the rate is reduced to values similar as the p ~ 
case above. 

Using the spin diffusion coefficient from Ref. [HJ for 
attractive interactions in the unitarity limit, we find good 
agreement with experiment both at low and intermedi- 
ate temperatures [lj. At high temperatures good agree- 
ment was also found r _1 ~ IODq/R 2 when the rate was 
reduced due to collisionless surface region [H|. At low 
temperatures the diffusion rate increases as T~ 2 but is 
reduced by the smaller inhomogeneity factor as observed 
in the experiments. The density dependence of the dif- 
fusion coefficient, which implies a difference in inhomo- 
geneity factor of ~ 6 as seen in Fig. (1), is therefore 
important in order to describe both the magnitude and 
temperature dependence of the diffusion rate. 

At the lowest temperatures accessible in the experi- 
ment of Sommer et al. superfluidity might start to 
appear in the center of the trap which is near spin equi- 
librium and therefore below the Clogston limit. As spin 
diffusion occurs in the center at low temperatures one 
may therefore look for increased spin diffusion and this 
way even determine the superfluid transition tempera- 
ture and the Clogston limit. 

The spin drag coefficient r s( j = ui 2 t measured in Ref. 
[l| can be expressed as 

hT sd = 2fiyTT? 

E F 3I p mD Q ' [ ' 
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It is independent of the trap frequencies as found in the 
experiments of Ref. The inhomogeneity factor affects 
the spin drag and the diffusion rate by increasing the ef- 
fective diffusion coefficient to D = I p Dq. The minimum 
of the diffusion coefficient Dq ~ h/m appears around 
medium temperature T ~ 0.5Tf and determines the uni- 
versal quantum limit to spin diffusivity in Fermi gases. 

In the unitarity limit at low temperatures the isother- 
mal conductivity k = 2Ep/3n(l + Fq) measured in Ref. 
[l| is compatible with the lower branch where the spin- 
symmetric Landau parameter Fq = /3 ~ —0.5. Likewise 
the spin susceptibility \ — ^Ep l^n{l+F§) is compatible 
with the lower branch spin-asymmetric Landau param- 
eter f a ~ 1 | Both are incompatible with the up- 
per branch where Fq is positive and Ffi negative. F.ex. 
in the Jastrow-Slater approximation [9] j3 — —0.54 in 
the unitarity limit on the lower branch but j3 — 2.93 
on the upper branch, and Fq — —7/3/5. Consequently 
the spin susceptibility is negative implying that the pre- 
ferred ground state is FM with phase separated spins, 
which contradicts the diffusion and mixing observed on 
longer time scales. 

In the experiments the two spin clouds are separated 
before ramping the magnetic field on resonance. When 
the barrier is removed and the two clouds come into con- 
tact the two-body wave function between unlike spins 
can evolve either into the lower or upper branches cor- 
responding to the BCS-BEC crossover or ferromagnetic 
limit respectively. On the longer diffusive time scales the 
experimental setup either favors the lower branch or the 
upper branch decays rapidly due to 3-body processes. If 
the FM state relaxes from the upper FM branch to the 
lower BCS branch before spin diffusion on the longer time 
scales, it would not necessarily contradict the claims in 
Ref. [10] that the two spin clouds are in a FM state 
at early times when they collide, bounce and undergo 
damped oscillations. 

In the recent work of Goulko et. al. [16j the Boltz- 
mann equation is solved at high temperatures off the 
unitarity limit. The bounce and oscillations of the two 
initially separated spin clouds and the initial dipole and 
breathing modes are calculated in detail. Viscous and 
diffusive relaxation are automatically included but mean 



fields are ignored in the high temperature approximation. 
Therefore the bounce and damped oscillations depend on 
the scattering cross section only which is symmetric for 
positive and negative scattering lengths. At lower tem- 
peratures we expect that the difference between the two 
scenarios of attractive or repulsive interactions becomes 
important, i.e. whether the two spin clouds evolve into 
the lower or upper branch with attractive or repulsive 
mean fields respectively. The collisional cross section 
may be mainly responsible for the bounce and oscilla- 
tion frequency but the sign should show up in the damp- 
ing. However, initially the two colliding spin clouds are 
far from mechanical and chemical equilibrium and the 
subsequent damped oscillations are dominated by ther- 
malization processes. The damped oscillations seem not 
to discriminate between the evolution towards a repul- 
sive or attractive initial state in the work by Goulko et 
al. as they find good agreement without attractive or re- 
pulsive mean fields. However, the calculations of Taylor 
et al. |l0] favor repulsive mean fields. Detailed analy- 
ses of the bounce, damped oscillations and diffusion are 
required at low temperatures before definite conclusions 
can be given. The colliding clouds are initially far from 
equilibrium and the thermalization affects the damped 
oscillations and diffusion. Especially a more accurate de- 
termination of the size of the overlap zone and its temper- 
ature would be useful for a better understanding of the 
damped oscillations and diffusion, since their rates are 
very temperature, density and spin density dependent. 

In summary, the strong density and temperature de- 
pendence of the spin diffusion coefficient leads to inho- 
mogeneous spin diffusion that takes place in the center 
of the cloud at low temperatures. This is opposite to the 
high temperature case where diffusion takes place in the 
surface and tail of the cloud circulating around the non- 
diffusive core. The calculated inhomogeneity factor for 
resonant attractive interactions leads to good agreement 
with experiments at low and higher temperatures but ex- 
cludes a FM phase at the longer diffusive time scales. The 
initial bounce of the two colliding spin clouds is far from 
mechanical and chemical equilibrium and the subsequent 
damped oscillations do not discriminate between the evo- 
lution towards a repulsive or attractive initial state. 
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